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ABSTRACT 


Simultaneous  measurements  of  waves,  and  vertical  and 
horizontal  water  particle  velocities  were  made  at  the 
breaker-line  within  the  surf  zone  using  a  capacitance  type 
penetrating  wave  staff,  a  pressure  wave  gauge,  and  an 
electromagnetic  current  meter.   Wave-  measurements  were  also 
made  at  seaward  and  shoreward  locations.   The  wave  energy- 
density  spectral  components  were  converted  to  velocity 
spectral  components  using  linear  wave  theory.   These 
computed  values  compared  well  qualitatively  with  the  meas- 
ured velocity  spectra.   Quantitatively,  the  results  showed 
that  linear  theory  underpredicted  wave-induced  horizontal 
velocity  spectral  components  by  about  50  percent  at  the 
frequency  of  peak  energy.   The  coherence  values  between 
waves  and  horizontal  velocity  were  high,  ranging  above  0.75. 
The  phase  angle  computation  showed  the  calculated  velocity 
components  leading  the  measured  velocity  components  by  an 
average  of  20  degrees,  indicating  an  unstable  wave  crest 
leading  the  particle  motion  in  the  body  of  the  wave. 
Probability  density  functions  were  computed  and  compared  to 
Gaussian  and  Gram-Charlier  distributions  using  the  chi-square 
goodness-of-fit  test.   The  Gram-Charlier  distribution 
qualitatively  gave  the  better  fit  to  the  data. 
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I.   INTRODUCTION 

A.   NATURE  OF  THE  PROBLEM 

The  problem  is  to  describe  the  kinematics  of  breaking 
waves  within  the  surf  zone  in  terms  of  the  water  surface 
elevation.   The  theories  which  have  been  developed  for  deep 
water  waves  can  be  applied  with  some  degree  of  certainty, 
and  can  be  tested  in  a  laboratory  situation.   These  theories 
generally  hold  up  to  the  time  that  the  wave  begins  to  shoal, 
and  with  somewhat  lesser  accuracy,  throughout  the  shoaling 
process  up  to  the  point  of  breaking.   At  the  breaker  point, 
however,  the  prediction  of  wave  kinematics  becomes  more 
difficult.   The  type  of  breaking  wave,  plunging,  surging,  or 
spilling,  depends  upon  deep  water  characteristics,  bottom 
type,  and  bottom  slope,  hence  requiring  a  different  descrip- 
tion of  the  kinematics  within  the  wave  for  each  breaker 
type.   Physical  limitations  preclude  the  construction  of  a 
realistic  laboratory  model  of  the  surf  zone.   And  finally, 
the  hostile  character  of  the  surf  zone  complicates  the 
taking  of  in  situ  measurements  to  the  degree  that  relatively 
little  data  exist  from  within  the  surf  zone. 

The  most  forthright  approach  to  the  problem  of  describing 
the  kinematics  of  breaking  waves  in  the  surf  zone  seems  to 
be  through  direct  measurements.   Advances  in  instrument  de- 
sign have  led  to  simple,  sturdy,  measuring  devices  with 
rapid  response  time  which  can  sense  the  small  scale  as  well 
as  the  large  scale  features.   By  making  the  assumption  that 
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the  processes  are  stationary,  this  assumption  being  valid 
at  least  over  a  short  record  of  20  to  30  minutes,  the  waves 
and  wave  induced  particle  velocities  can  be  dealt  with  as 
stationary  random  phenomena,  thereby  allowing  the  utiliza- 
tion of  standard  statistical  analysis  procedures. 

B.   REVIEW  OF  PREVIOUS  WORKS 

The  study  of  the  kinematics  of  breaking  waves  in  the 
surf  zone  has  progressed  slowly  both  due  to  the  difficulty 
in  making  direct  field  measurements,  and  due  to  the  diffi- 
culty in  modeling  the  surf  zone  in  the  laboratory. 

Early  work  by  Iversen  (1953)  involved  use  of  a  wave 
channel  in  which  waves  were  generated  using  a  hinged  plane 
oscillating  flap.  Measurements  were  made  using  photographic 
techniques.  The  length  of  the  channel  limited  the  wave  type 
to  plunging  and  surging  breakers.  Adeyemo  (1970)  made  simi- 
lar laboratory  measurements  using  hydrogen  bubble  and  photo- 
graphic methods. 

Inman  (1956)  made  field  measurements  of  water  particle 
velocity  by  measuring  the  drag  force  under  the  wave  in  order 
to  infer  particle  velocities.   Miller  and  Zeigler  (1964)  used 
both  acoustic  and  electromagnetic  current  meters  to  determine 
the  particle  motion  in  the  surf  zone.  Walker  (1969)  made 
studies  using  propeller  type  flow  meters.   Thornton  (1968) 
and  Steer  (1972)  made  use  of  electromagnetic  current  meters 
to  measure  the  particle  motion. 


C.   OBJECTIVE 

The  objective  of  this  research  is  to  study  the  kinematics 
of  water  particles  in  breaking  waves  within  the  surf  zone. 
Simultaneous  measurements  were  made  of  the  instantaneous 
water  surface  elevation  and  two  orthogonal  water  particle 
velocities  at  a  fixed  location  in  the  surf  zone.   Additional 
wave  measurements  were  made  at  points  seaward  and  shoreward 
of  the  primary  location.   Estimates  of  the  probability 
density  functions  and  power  spectra  of  the  wave  height  and 
particle  velocities  are  to  be  made.   The  applicability  of 
linear  wave  theory  is  to  be  measured  in  computing  velocity 
spectral  components  using  the  energy-density  spectrum  of 
the  waves.   The  computed  velocity  spectra  are  to  be  compared 
with  actually  measured  velocity  spectra. 
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II.   STATISTICAL  ANALYSIS 

It  was  assumed  that  the  data  to  be  analyzed  represents 
realizations  of  stationary,  random  processes.   This  assump- 
tion allows  the  use  of  standard  statistical  techniques.   The 
basic  analyses  are  described  in  the  following  sections  as 
if  performed  on  continuous  data  of  infinite  duration  for 
ease  of  discussion.   The  actual  analyses  were  performed 
digitally  on  records  of  finite  length. 

A.   PROBABILITY  DENSITY  FUNCTION 

The  probability  density  function  describes  the  probability 
that  a  sample  time  history  record  x(t)  will  assume  a  value 
between  x  and  (x  +  Ax) ,  and  that  this  approaches  an  exact 
probability  as  T,  the  observation  time,  approaches  infinity, 
and  as  the  defined  range,  Ax,  approaches  zero.   The  proba- 
bility density  function  is  given  mathematically  as: 


T 
P(x)  =  Lim  Lim  ~-  -±  (2.1) 

Ax+0  T+~   ax   x 


Where  T.  is  the  amount  of  time  that  x(t)  falls  within  the 
1 

range  (x,  x  +  Ax) . 

The  computed  probability  density  function  was  compared 
to  the  Gaussian  and  Gram-Charlier  distributions  and  tested 
for  comparison  using  the  chi-square  goodness-of-fit  test. 
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The  Gaussian  probability  density  function  is  given  as: 


P   (X)  -  1 e-(x-y)2/2a2  (?2) 

Ga      c  V2T 


where   a  =  standard  deviation 
y  =  mean 

and  in  standardized  form,  with  z  =  ^—^-   , 

P  (z)  =_L_e-2/2 
Ga     VTfT 

The  Gaussian  distribution  is  completely  described  by  the 
mean  and  variance,  and  has  zero  skewness  and  kurtosis  equal 
to  three . 

The  Gram-Charlier  distribution  was  shown  by  Longuet- 
Higgins  (1963)  to  better  describe  the  probability  distribu- 
tion of  the  instantaneous  water  surface  elevation,  f\,    in 
moderately  non-linear  swell  in  deep  water.   It  is  given  by: 

m  m^  m 

PGc(z)    =    PGa(z)[1'°   +    2?  H3   +    (27H4    +  TT  V    +    •••]         (2"3) 


where 


H     -   „n       n(n-l)    z11"2    ,    n(n-l)  (n-2)  (n-3)    zn~4 
n        z      "        1!  2  2!  22 

(Hermite  Polynomial) 

m~  =  skewness  of  record 

nu  =  kurtosis  -  3 
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It  can  be  seen  that  for  zero  skewness  and  kurtosis  of  three, 
the  fourth  order  Gram-Charlier  distribution  degenerates  to 
the  Gaussian  distribution. 

To  apply  the  chi-square  goodness-of-fit  test,  the  random 
variable,  x(t),  was  divided  into  K  intervals,  called  class 
intervals,  thereby  forming  a  frequency  histogram.   The 
observed  frequency,  f . ,  is  the  number  of  observations  within 
the  ith  class  interval.   The  expected  frequency,  F. ,  is  the 
number  of  observations  which  fall  in  the  ith  interval  of  the 
theoretical  distribution.   The  total  discrepancy  between  the 
observed  and  expected  frequency  is  given  by: 

„    K   (f.  -F.  )2 
X2  =  Z  \      X  (2.4) 

i=l     1 

2 

As  the  value  of  x   increases  from  zero,  the  goodness-of- 
fit  decreases. 

B.   SPECTRAL  ANALYSIS 

1 .   Energy-Density  Spectrum 

The  auto-covariance  for  a  record,  x(t),  is  defined 
mathematically  by 


,    t+T 
R   (t)  =  Lim  ±   /    x  (t)x  (t+T)dt         (2.5) 
x±  t->°°     t 


The  value  obtained  describes  the  dependence  of  the  values  of 
the  data  set  at  one  point  in  time  on  the  values  at  some 
later  time,  f,  inthe  dataset.  The  function  is  real-valued  and 


13 


even,  and  has  a  maximum  at  t  =  0.   The  above  expression 
implies  an  infinitely  long  data  set.   To  deal  with  the 
reality  of  a  data  set  of  finite  length,  a  lag  window  is 
applied.   The  choice  of  lag  window  for  this  analysis  was  the 
Parzen  lag  window,  because  it  has  the  advantage  of  no 
negative  side  lobes,  which  maintains  the  coherence  function 
between  its  theoretical  values  of  +  1.   It  also  insures 
no  numerical  resonance  problems  when  computing  cross-spectra. 
The  Parzen  window  is  given  by  Bendat  and  Piersol  (1966)  as: 


T 

1  -  6(^)2  +  6  (^)3   ,  x  =  0,1,2,...,^ 

m         m 
P(t)  t  o  T  (2.6) 

2  M  -  — V  t  =  —  +  1      T 

^   L-L      rp   J  ,  T       ^   T   ±,  .  .  .  ,  lm 

m 

}       m 


The  energy-density  spectrum  is  computed  from  the  Fourier 
transform  of  the  covariance  function  as  modified  by  the 
Parzen  lag  window.   The  spectrum  is  given  by: 


S11(f)  =   /   P(T)R11(T)e-i27TfTdT  (2.7) 


where  f  is  a  particular  frequency  of  interest. 
2 .   Cross-Spectral  Density 

The  cross-covariance  function  describes  the  general 
dependence  of  the  values  of  one  set  of  data  on  the  other  set, 
in  contrast  to  the  auto-covariance  function  which  compares  a 
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record  against  itself.   The  cross-covariance  function  is 
given  mathematically  as: 

..   t+T 
R   (t)  =  Lim  -   /    x  (t)x?(t+x)dt  (2.8) 

The  quantity  R, p(x)  is  a  real-valued  function,  but  it  is  not 
an  even  function,  and  does  not  necessarily  have  a  maximum 
value  at  x  =  0. 

A  lag  window  again  is  applied  to  the  cross- 
covariance  function  before  computing  the  cross-spectral 
density.   Applying  the  Parzen  window  and  taking  the  Fourier 
transform,  the  cross-spectral  density  is  given  by: 


S   (f)  =   /  P(x)R12(T)e"i2Tr;fTdT  (2.9) 

—  00 


Since,  unlike  the  energy-density  spectrum,  the  cross-spectrum 
is  not  an  even  function,  it  is  expressed  in  real  and 
imaginary  parts  by: 


S12(f)  =  C12(f)  -  i  Q12(f)  (2.10) 


The  real  part,  called  the  co-spectrum,  is  given  by: 


C12(f)  =  2   /  P(t)[R12(t)  +  R12(-T)]cos(2iTfT)dT  (2.11) 
The  imaginary  part,  called  the  quadrature  spectrum,  is  given 


by: 


Q12(f)  =  2   /  P(t)[R12(t)  +  Ri2(-T)]sin(27TfT)dT  (2.12) 
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3.   Phase  Angle 

Equation  (2.10)  can  be  rewritten  in  complex  polar 
notation  in  terms  of  the  modulus  of  the  cross-spectral 
density  and  a  phase  angle, 

S12(f)  =  |S12(f)|e~ie12(f)  (2.13) 

where  e,p(f)  =  average  angular  difference  by  which  the 
cross-correlated  components  of  x?(t)  lead  those  of  x, (t)  in 
each  spectral  band.   The  phase  angle  expressed  in  terms  of  : 
the  co-spectra  and  quadrature-spectra  is  given  by: 

e19(f)  =  tan  1[-^ ]  (2.14) 

12  C   (f) 


where 


-it  <  eno(f)  <  7T 


.c. 


4.   Coherence 

An  expression  for  the  comparison  of  two  random 
variables  as  a  function  of  frequency  is  described  by  the 
coherency  squared,  or  simply  coherence,  and  is  derived  using 
the  expressions  previously  stated  for  energy-density  spectra 
and  for  cross-spectral  density.   In  mathematical  terms,  the 
coherence  function  is  given  as: 

|s12(f)|2 

Y-,?2(f)  =  (2.15) 

ld  su(f)s22(f) 


where : 


0  <  ^122  <  1 
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The  coherence  function  has  the  characteristics  that  when 

2 

Y-,  p  (f)  =  1,  the  two  records  are  completely  correlated  at 

2 
that  frequency;  and  when  y1?  (f)  =  0 ,  the  records  are 

statistically  independent  at  that  frequency. 


17 


III.   EXPERIMENTAL  PROCEDURE 

A.  EXPERIMENTAL  SITE 

The  experiment  was  conducted  at  Del  Monte  Beach  inside 
Monterey  Bay  adjacent  to  Monterey,  California,  on  12  April, 
1973.   The  beach  profile  shown  in  Figure  (1)  is  one  of  gentle 
slope.   Plunging  and  spilling  breakers  are  very  common  at 
this  site,  with  spilling  breakers  being  the  predominant  type. 
The  waves  at  this  location  are  usually  severely  directionally 
filtered  by  refraction  due  to  the  geometry  of  the  coastline 
and  bay,  and  approach  almost  perpendicular  to  the  shore. 

B.  INSTRUMENTATION 

The  primary  measuring  instruments  were  placed  on  a  three 
meter  steel  tower  as  shown  in  Figure  (2).   The  tower  was 
guyed  so  that  tower  movement  and  vibration  were  negligible. 
The  instruments  were  arranged  so  as  to  be  in  the  same  verti- 
cal plane,  with  the  exception  of  a  small  offset  parallel  to 
the  wave  crests  given  to  the  water  current  meter  to  prevent 
interaction  with  the  capacitance  probe. 

The  water  particle  velocities  were  measured  with  an 
Engineering  Physics  Company  water  current  meter  Model  6130. 
The  sea  surface  elevation  was  measured  with  a  capacitance 
type  penetrating  wave  staff  fabricated  locally,  and  with 
an  Interstate  Electronics  Model  DP-200  portable  wave  recor- 
der with  SDP-201  differential  pressure  sensor. 
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Additional  instrumentation  was  placed  seaward  and  shore- 
ward of  the  main  tower.  A  pressure  wave  gauge  was  installed 
seaward  of  the  main  tower  by  mounting  it  on  a  weighted  sled. 
The  sled  was  towed  into  position  by  attaching  a  winch  cable 
from  R.  V.  Acania.  Acania,  in  her  anchored  position,  heaved 
around,  pulling  the  sled  into  position.  The  mean  water  at 
the  sled  during  the  experiment  was  2.69  meters. 

A  Baylor  Company  Wave  Staff  System  Model-  13528R  was 
installed  12.9  meters  shoreward  on  a  steel  tower.   This  wave 
staff  was  placed  in  a  line  with  the  primary  tower  and  the 
sled  so  as  to  be  on  an  orthogonal  with  the  predicted  wave 
crests . 

The  capacitance  wave  staff  was  designed  by  L.  F. 
McGoldrick  (1969)  at  the  University  of  Chicago  as  an  inex- 
pensive, but  ultra-linear  highly  sensitive  instrument  for 
the  detection  of  capillary  waves  in  the  20  to  30  HZ  range. 
It  operates  on  the  principal  that  changing  the  plate  dimen- 
sions of  a  capacitor  changes  its  capacitance  and  hence  its 
influence  on  an  electrical  circuit.   The  probe  was  constructed 
of  one-eighth  inch  0.  D.  bronze  wire;  a  dialectric  of  one- 
sixteenth  inch  wall  thickness  polystyrene  tubing  which 
fitted  tightly  over  the  wire;  and  appropriate  electrical 
connections.   The  sea  water  in  which  the  probe  was  immersed 
constituted  the  second  plate  of  the  capacitor.   The  probe 
was  held  stationary  on  the  tower,  and  as  the  sea  level  fluc- 
tuated, the  capacitance  changed  accordingly.   These  fluctua- 
tions were  sensed  by  a  transistorized  circuit  operated  by 
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two  self-contained  ten  volt  batteries.   The  electronics 
section  was  housed  in  a  watertight  case  which  was  mounted 
on  the  tower.   This  enabled  the  connecting  leads  to  be 
relatively  short,  thereby  minimizing  wire-to-wire  capacitance 
A  static  calibration  of  the  system  was  made  by  immersing 
the  staff  in  a  water  tank.   The  resulting  calibration  is 
shown  in  Figure  (3). 

The  Interstate  Electronics  Pressure  Sensor  used  on  the 
primary  tower  and  at  the  sled  site  is  a  small,  unbonded 
strain  gauge  bridge.   The  sea  pressure  is  coupled  by  a  neo- 
prene  diaphragm  to  a  silicone  fluid  filling  the  interior. 
The  arrangement  is  such  that  slow  changes,  such  as  tides, 
are  lost  through  the  hydraulic  filter,  and  the  transducer 
only  senses  the  rapid  pressure  fluctuations  of  waves.   The 
high-pass  filter  response  time  constant  is  ^8  seconds. 

The  Engineering  Physics  Company  electromagnetic  current 
meter  measures  two  orthogonal  components  of  water  particle 
velocity  through  a  range  of  zero  to  one  and  one-half  M/SEC 
with  a  maximum  output  error  of  one  percent  of  full  scale 
reading.   The  instrument  has  two  electrical  time  constant 
settings,  0.3  and  1.0  seconds.   The  operation  of  the  unit 
is  based  on  Faraday's  principle  of  magnetic  induction. 
Faraday  discovered  that  currents  are  induced  in  a  conductor 
which  is  put  in  motion  through  a  magnetic  field.   This 
induction  can  be  expressed  by  the  vector  equation: 

E  =  U  X  B 
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where  e  is  the  induced  electric  field,  U  is  the  velocity  of 

-*- 
motion,  and  B  is  the  intensity  of  the  magnetic  induction. 

In  the  sea  the  induced  electric  field,  E,  is  linearly  pro- 
portional to  the  amount  of  fluid  flowing  through  the  magnetic 
field.   The  magnetic  field,  B,  is  most  intense  near  the  probe, 
and  decreases  with  the  distance  from  the  probe  according 
to  the  inverse  square  law.   Thus,  the  flow  velocity  near  the 
probe  is  weighted  more  heavily  than  that  away  from  the  probe, 
with  only  that  water  which  is  within  a  distance  of  two  to 
three  probe  radii  contributing  significantly  to  the  flow 
induced  voltage.   The  radius  of  the  probe  is  3/8  inch. 

The  flow  meter  was  recalibrated  using  the  method  of  Kra- 
pohl  (1972)  and  Steer  (1972).   This  procedure  involves  oscil- 
lating the  probe  in  a  water  tank.   The  equipment  arrangement 
is  shown  in  Figure  (4).   The  flow  meter  probe  was  mounted 
on  a  carriage  which  travelled  back  and  forth  on  rails.   The 
carriage  was  driven  by  a  variable  speed  motor  geared  to  an 
eccentric  throw  arm.   The  peak  carriage  velocity  was  calcu- 
lated to  be  the  tangential  velocity  of  the  throw  arm.   The 
values  of  maximum  carriage  velocity,  V, ,  and  maximum  measured 
velocity,  V  ,  were  determined  for  different  angular  frequen- 
cies.  Each  pair  of  electrodes  was  calibrated  by  orienting 
them  parallel  to  the  flow. 

Figure  (5)  is  a  plot  of  the  ratio  of  measured  maximum 
velocity  and  maximum  carriage  velocity  vs.  frequency.   The 
instrument's  response  is  shown  to  be  flat  out  to  frequencies 
of  about  0.5  Hz.   At  higher  frequencies,  the  cart  became 
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unstable.   In  addition,  waves  were  generated  in  the  tank 
which  affected  the  flow  pattern.   This  upper  frequency  limi- 
tation was  net  deemed  a  significant  factor  in  measurement 
of  ocean  waves  in  which  the  significant  energy  was  found  to 
be  below  0.5  Hz. 

The  Model  13529R  Baylor  Company  Wave  Staff  System  is  an 
instantaneous  water  level  measuring  device  consisting  of  two 
tensional  half  inch,  6x19  IWRC  SS  wire  ropes.   Together  with 
these  lengths  of  stainless  steel  wire  rope,  the  transducer 
produces  an  electrically  linear  direct  current  output  that 
is  proportional  to  the  amount  of  wave  staff  above  a  short 
circuit  produced  by  the  water  surface.   The  working  length 
of  the  staff,  through  previous  modifications,  had  been 
shortened  from  50  to  5  feet.   The  results  of  the  recalibra- 
tion  are  shown  in  Figure  (6). 

C.   DATA  COLLECTION 

The  experiment  was  performed  on  12  April,  1973.   The 
towers  and  sled  had  been  installed  at  low  tide  two  days  pre- 
viously.  Tests  were  conducted  and  final  calibration  checks 
were  made  on  11  and  12  April.   Included  in  the  calibration 
was  the  placing  of  a  two  volt  D.  C.  reference  voltage  on 
the  FM  tape  to  be  used  during  data  processing.   This  insured 
that  the  signal  gains  could  be  checked  after  digitizing. 

During  the  afternoon  of  12  April,  swell  from  a  storm  off 
the  Washington  coast  began  to  arrive,  with  a  buildup  to  peak 
wave  height  coincident  with  high  tide.   During  the  experimer  , 
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the  significant  wave  height  at  breaking  was  about  1  meter. 
The  waves  were  generally  of  the  plunging  type,  and  were 
breaking  in  the  general  area  of  the  main  tower.   The  breaker 
point  was  visually  observed  to  move  seaward  and  shoreward 
according  to  the  incident  wave  height. 

All  data  were  simultaneously  recorded  using  a  Sangamo 
Model  3500 j  1^  channel  FM  magnetic  tape  recorder.   Continuous 
monitoring  of  data  was  accomplished  using  strip  chart 
recorders  and  voltmeters. 

D.   DATA  PRE-PROCESSING 

The  data  were  transcribed  using  the  CS-500  Analog  Compu- 
ter to  an  eight  channel  Clevit  Brush  strip  chart  recorder 
for  signal  verification.   It  was  digitized  and  recorded  onto 
seven  track  tape  using  the  XDS-9300  computer  in  conjunction 
with  the  CS-500  Analog  Computer.   The  digitized  records  were 
transferred  to  nine  track  tape  for  analysis  by  the  IBM  36O 
computer. 
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IV.   ANALYSIS  OF  DATA 

Figures  (7)  and  (8)  are  typical  records  showing  that 
the  breaking  waves  were  roughly  triangular  in  shape  and  were 
very  peaked.   This  corresponds  to  the  general  observation 
that  plunging  breakers  at  or  near  the  breaker  point  have  a 
steep,  almost  vertical  front  face,  and  a  sloping  profile 
toward  their  trailing  edge.   The  records  show  small  peaks 
between  the  primary  wave  peaks.   It  is  thought  that  the 
small  peaks  are  harmonics  of  the  primary  wave  frequency. 

The  capacitance  gauge  showed  a  more  sharply  peaked  wave 
than  did  the  pressure  gauges.   This  was  due  to  the  nature  of 
the  instruments  in  that  the  capacitance  gauge  could  sense 
more  rapid  changes  in  the  wave  slope. 

A.   SAMPLING  PROCEDURE 

The  data  were  read  at  a  rate  of  3  samples  per  second, 
corresponding  to  a  sample  interval  of  0.33  seconds.   This 
resulted  in  a  Nyquist  frequency  of  1.5  Hz.   The  sampling 
rate  was  considered  to  be  sufficiently  high  to  avoid  aliasing 
of  energy  into  the  portion  of  the  spectra  which  was  of 
interest  from  observations  of  the  analog  record. 

The  length  of  record  analyzed  was  21  minutes.   The  maxi- 
mum lag  time  was  taken  as  five  percent  of  the  record  giving 
a  spectral  bandwidth  resolution  of  O.OO873  Hz.   According  to 
Blackman  and  Tukey  (1958),  this  yields  40  degrees  of  freedom. 
By  applying  the  chi-square  distribution,  the  80  percent 
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confidence  limits  for  40  degrees  of  freedom  are  found  to 
be  between  0.77  and  I.38  of  the  measured  power  spectral 
estimates . 

B.   PROBABILITY  DENSITY  FUNCTION 

The  probability  density  function,  pdf,  for  each  wave 
and  water  particle  velocity  record  was  estimated  using  the 
procedures  outlined  in  Section  II-A.   For  each  record,  the 
standard  deviation,  variance,  skewness,  and  kurtosis  were 
computed  and  are  summarized  in  Table  I.   Measured  distribu- 
tions were  compared  with  the  Gaussian  and  Gram-Charlier 
distributions  using  the  chi-square  goodness-of-fit  test. 
The  various  measured  distributions  were  also  plotted  against 
these  analytical  distributions. 

1 .   Sea  Surface  Elevation 

At  the  seaward  station,  the  measured  wave  pdf  shows 
some  relation  to  the  Gaussian  distribution  [Figure (9)]. 
The  Gram-Charlier  gives  a  better  fit  in  terms  of  the  chi- 
square  goodness-of-fit  test.   This  is  to  be  expected  in 
that  the  Gram-Charlier  pdf  is  calculated  using  the  assigned 
parameters  from  the  data,  and  including  moments  higher  than 
the  second.  . 

At  the  primary  station,  the  pdf  again  follows  the 
Gram-Charlier  more  closely  than  the  Gaussian,  but  not  to  as 
great  a  degree  as  in  the  deeper  water  case  [Figures  (10) 
and  (11)]. 
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At  the  shoreward  station,  the  distribution  again 
more  closely  resembles  the  Gram-Charlier  [Figure  (12)]. 

The  positive  skewness  measured  at  each  station  indi- 
cates a  large  number  of  small  negative  values  and  a  higher 
probability  of  large  positive  values  of  wave  height  and 
velocity  when  compared  to  the  Gaussian.   This  qualitative 
observation  supports  the  notion  that  in  shallow  water,  the 
wave  crests  are  steep  and  narrow,  and  the  wave  troughs  are 
flat  and  wide . 

2  .   Water  Particle  Velocity 

The  probability  density  functions  shown  in  Figures 
(13)  and  (14)  for  water  particle  velocity  again  follow  the 
Gram-Charlier  distribution  more  closely  than  they  do  the 
Gaussian.   The  Gram-Charlier  distribution  was  computed  by 
carrying  the  series  to  include  the  fourth  moments  of  the 
frequency  distribution  (Equation  2.3).   This  appeared  to 
be  sufficient  at  the  seaward  station  [Figure  (9)].   However, 
due  to  the  more  non-linear  situation  in  the  surf  zone,  it 
appears  that  higher  moments  should  be  included.   This  would 
give  a  better  fit  and  get  rid  of  the  humps  in  the  tails  of 
the  distribution. 

C.   SPECTRAL  ANALYSIS 

The  energy-density  spectrum  was  calculated  for  each  of 
the  data  records.   Cross-spectra  were  then  computed  between 
variables,  from  which  the  coherence  and  phase  were  determined 
Wave  energy-density  spectral  components  were  computed  from 
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pressure  energy-density  components  using  linear  wave  theory. 
Linear  wave  theory  was  also  used  to  transform  the  wave  spec- 
tral components  to  spectral  estimates  of  horizontal  and  ver- 
tical velocities.   These  calculated  spectral  values  were 
then  compared  to  measured  particle  velocity  spectra. 
1.   V/ave  Height 

The  capacitance  wave  gauge  was  taken  as  the  primary 
wave  height  measuring  instrument.   The  capacitance  gauge  gave 
a  direct  measurement  of  instantaneous  sea  surface  elevation, 
and  by  taking  the  mean,  a  measure  of  the  average  sea  surface 
elevation  was  obtained. 

The  waves  were  also  measured  using  a  pressure  wave 
gauge.   To  spectrally  convert  these  pressure  fluctuations  to 
values  of  water  surface  elevation,  a  transfer  function  was 
applied.   This  function  is  based  on  linear  theory,  and  is 
given  as: 

cosh  kh     2  (/j,1) 

S(f)(wave)  =  ^"oshklh+z)-1    ^  ']  (pressure  ) 

where:   S(f)  =  energy-density  spectra 

k  =  wave  number 

h  =  water  depth 

z  =  pressure  sensor  depth 

The  energy-density  spectra  of  both  the  capacitance 
and  pressure  wave  gauges  are  shown  in  Figure  (15).   The 
pressure  spectrum  has  been  converted  to  a  wave  spectrum 
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using  Equation  U . 1 .   The  spectra,  which  are  characteristic 
of  all  the  records,  show  a  narrow  banded  energy-density  peak 
at  a  frequency  of  0.07  Hz.,  corresponding  to  a  period  of  l!| .  3 
seconds.   Also  evident  are  peaks  at  0.1'J  and  0.21  Hz.  which 
appear  to  be  harmonics  of  the  primary  peak  at  0.07  Hz.   These 
harmonics  appear  to  be  physical  as  was  observed  from  the  strip 
chart  earlier.   These  harmonics  probably  also  have  energy 
contributions  due  to  the  Fourier  computational  technique.   A 
low  frequency  peak  is  also  shown  at  about  0.02  Hz.  corres- 
ponding to  a  period  of  50  seconds.   Energy  at  this  frequency 
is  commonly  measured  along  this  coast  and  is  attributed  to 
surf  beat. 

A  peak  occurs  in  all  spectra  at  0.6  Hz.   This  is 
attributed  to  mechanical  noise  in  the  tape  recorder.   This 
peak  is  at  a  frequency  outside  the  range  of  appreciable  wave 
energy  and  does  not  affect  the  results  of  the  experiment. 

The  magnitude  of  the  wave  spectral  components  at 
the  frequency  of  peak  energy  density  calculated  from  the 
pressure  record  using  linear  wave  theory  was  approximately 
60  percent  of  that  measured  by  the  capacitance  wave  gauge. 
The  coherence  values  were  very  high,  ranging  to  greater 
then  0.9  in  the  maximum  energy  portion  of  the  wave  spectrum. 
The  phase  spectrum  shows  the  capacitance  wave  gauge  measure- 
ments generally  leading  the  pressure  measurements  by  a  few 
degrees,  probably  due  to  the  ability  of  the  capacitance 
gauge  to  better  measure  the  curling  portion  of  the  wave. 


M 


The  spectra  of  the  offshore  wave  gauge  and  the  capac- 
itance wave  gauge  were  compared,  as  were  the  inshore  and 
capacitance  gauges  [Figures  (16)  and  (17)].   A  dynamic  pressure 
correction  [Equation  4.1]  was  not  applied  to  the  seaward 
pressure  gauge.   The  results  show  that  the  waves  were 
shoaling  throughout  the  process;  that  is,  there  was  a  gradu- 
al buildup  of  energy  as  the  wave  progressed  from  outside 
the  breaker  zone  into  the  breaker  zone . 

The  greater  variance  at  the  inshore  wave  gauge  com- 
pared with  the  wave  gauge  at  the  breaker  line  is  attributed 
to  energy  reflection  from  the  beach.   The  Baylor  wave  gauge 
tower  was  approximately  25  meters  from  the  point  at  which 
mean  water  level  intersects  with  the  beach,  and  12.5  meters 
shoreward  of  the  mean  breaker  line  where  the  flow  meter  and 
primary  wave  gauges  were  located.   These  distances  correspond 
to  one-half  and  one-quarter  wave  lengths,  respectively,  of 
a  13  second  period  progressive  wave  in  a  depth  of  1.9  meters. 
This  arrangement  places  a  node  at  the  breaker  line  and  an 
anti-node  at  the  Baylor  wave  gauge  tower  for  a  standing  wave 
of  this  period.   This  13  second  period  corresponds  closely 
to  the  period  of  peak  energy. 

The  phase  for  a  theoretical  progressive  wave  is 
plotted  with  the  measured  phase  for  comparison  in  Figure  (17). 
The  phase  relationship  was  calculated  using  linear  wave 
theory.   Deviations  from  the  calculated  progressive  phase 
angle  could  be  attributed  to  wave  reflections. 
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2 .   Sea  Surface  Elevation  and  Water  Particle  Velocity 
The  spectra  of  sea  surface  elevation  and  water 
particle  velocities  were  compared  by  using  linear  wave 
theory  to  convert  from  wave  spectral  components  to  water 
particle  velocity  spectral  components.   The  water  surface 
elevation  described  by  linear  theory  is  given  by: 

n  =  a  cos  (kx-at)  (4.2) 


where 


n  =  sea  surface  elevation 

a  =  wave  amplitude 

k  =  wave  number 

x  =  horizontal  cartesian  coordinate 

a   =   frequency 

t  =  time 


The  velocity  components  are  given  by: 

n      cosh  k(h  +  z)  nr.^(VY  -  est)    (U    1) 
Horizontal  velocity  =  u  =  aa  — slnh  kh cosux  at;  ^.j; 


«   slnh  k(h  +  z)   ,  (v     _     t \    ru    jn 
Vertical  velocity   =  w  =  ao  — ^nTTkh    smux  oz)    K*-*) 


The  horizontal  velocity  spectrum  calculated  from  the  wave 
spectrum  is  given  by: 


o/*.n        r   cosh  k(h  +  z)n2  ~(f\  (i|.5) 

S(f)(vel)  =  [a  "iinlTlS    J   S(f)(wave) 


and  the  vertical  velocity  spectrum  by 
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qffM       -r«sinhk(h  +  z)-,2  ,    - 

S(f)(vel)  -  C°  -iTnh  kh   ]   S(f>(wave)  ^-6) 

where  the  terms  in  brackets  are  the  corresponding  transfer 
functions . 

a.  Horizontal  V/ater  Particle  Velocity 

Linear  theory  under-calculates  the  horizontal 
velocity  spectrum  by  50  percent  at  the  peak  energy  point 
[See  Figure  (18)].   The  values  of  coherence  are  high, 
ranging  above  0.75.   The  phase  angle,  which  according  to 
linear  theory  should  equal  zero,  ranges  up  to  4 3  degrees 
from  its  average  value  of  about  20  degrees.   This  shows  the 
crest  of  the  unstable  breaking  wave  leading  the  horizontal 
particle  velocities  in  the  body  of  the  flow.   This  compares 
well  with  the  laboratory  findings  by  Adeyemo  (19  70)  that 
the  maximum  shoreward  horizontal  velocity  for  monochromatic 
waves  occurs  25  degrees  after  the  passage  of  the  wave  crest. 

b.  Vertical  Water  Particle  Velocity 

The  measured  rms  vertical  velocity  was  an  order 
of  magnitude  less  than  the  measured  rms  horizontal  velocity, 
and  showed  very  little  wave  induced  motion  [Figure  (19)]. 
The  coherence  reached  a  maximum  value  of  0.37  at  a  frequency 
of  0.15  Hz.,  which  is  above  the  frequency  of  maximum  energy 
as  measured  for  the  wave  and  horizontal  velocity  spectra. 
A  phase  angle  of  130°  to  145°  was  computed,  in  contrast  to 
the  theoretical  value  of  90°  between  wave  height  and  vertical 
velocity . 
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The  lack  of  coherence  between  water  surface 
elevation  and  vertical  velocity  could  be  due  to  a  noisy 
electrical  signal  from  the  instrument,  a  conversion  of 
wave-induced  motion  to  turbulent  motion,  or  both. 
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V.   CONCLUSIONS 

The  probability  density  functions  for  the  various  records 
compared  better  with  the  Gram-Charlier  distribution  than  with 
the  Gaussian  distribution  when  tested  using  the  chi-square 
goodness-of-fit  test.   This  result  verifies  the  observed 
asymmetrical  shape  of  the  waves  in  the  breaker  zone  and 
points  out  the  importance  of  including  higher  order  moments 
in  describing  wave  phenomena  in  very  shallow  water. 

The  values  of  the  horizontal  velocity  energy-density 
spectral  components  calculated  from  wave  spectra  using  linear 
theory  indicate  a  qualitative,  but  not  a  quantitative  rela- 
tionship.  In  general,  linear  theory  underpredicts  the 
magnitude  of  the  velocity  spectra. 

The  high  coherence  between  waves  and  horizontal  velocity 
shows  that  most  of  the  motion  in  the  body  of  the  breaking 
wave  is  wave-induced.   A  similar  conclusion  could  not  be 
drawn  when  comparing  waves  and  vertical  velocity. 

The  spectral  difference  between  water  surface  elevation 
and  horizontal  velocity  was  near  that  predicted  by  linear 
theory.   The  waves  lead  the  horizontal  velocity  in  spectral 
phase  on  the  average  by  20  degrees,  implying  that  the 
"curling"  crest  of  the  wave  arrives  prior  to  maximum  water 
particle  velocity. 

In  summary,  the  water  particle  motion  in  the  breaking 
waves  seems  to  be  very  much  wave  induced  in  the  body  of  the 


53 


flow.   Linear  theory  under-predicts  the  magnitude  of  the 
velocity  spectra  when  computing  from  wave  spectra. 
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